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ABSTRACT 

We present the results of a systematic Green Bank Telescope (GBT) and Giant Metre- 
wave Radio Telescope (GMRT) survey for 21-cm absorption in a sample of 10 damped 
Lyman-a; (DLA) systems at 2 < z a bs< 3.4. Analysis of L-band Very Long Baseline 
Array (VLBA) images of the background QSOs are also presented. We detect 21-cm 
absorption in only one DLA (at z a b s = 3.1745 towards J1337+3152). Thus the detec- 
tion rate of 21-cm absorption is ~ 10% when no limit on the integrated optical depth 
(/ r(v)dv) is imposed and ~ 13% for a 3a limit of 0.4 km s _1 . Combining our data 
with the data from the literature (a sample of 28 DLAs) and assuming the measured 
core fraction at milliarcsecond scale to represent the gas covering factor, we find that 
the H 1 gas in DLAs at z > 2 is predominantly constituted by warm neutral medium. 
The detection rate of 21-cm absorption seems to be higher for systems with higher 
N(H 1) or metallicity. However, no clear correlation is found between the integrated 
21-cm optical depth (or the spin-temperature, Tg) and either N(H 1), metallicity or 
velocity spread of the low ionization species. There are 13 DLAs in our sample for 
which high resolution optical spectra covering the expected wavelength range of H2 
absorption are available. We report the detection of H2 molecules in the z a b s = 3.3871 
21-cm absorber towards J0203+1134 (PKS 0201+113). In 8 cases, neither H 2 (with 
molecular fraction f(H.2)< 10~ 6 ) nor 21-cm absorption (with Tg// C > 700 K) are 
detected. The lack of 21-cm and H2 absorption in these systems can be explained if 
most of the H 1 in these DLAs originate from low density high temperature gas. In one 
case we have a DLA with 21-cm absorption not showing H2 absorption. In two cases, 
both species are detected but do not originate from the same velocity component. In 
the remaining 2 cases 21-cm absorption is not detected despite the presence of H2 
with evidence for the presence of cold gas. All this is consistent with the idea that the 
H2 components seen in DLAs are compact (with sizes of < 15 pc) and contain only a 
small fraction (i.e typically < 10%) of the total N(H 1) measured in the DLAs. This 
implies that the molecular fractions /(H2) reported from the H2 surveys should be 
considered as conservative lower limits for the H2 components. 
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The Galactic interstellar medium (ISM) has a multiphase 
structure with neutral hydrogen being distributed between 
the cold neutral (CNM), warm neutral (WNM) and warm 
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1 INTRODUCTION 



ionized (WIM) media. A large fraction of the gas is also 
found in diffuse, translucent and dense molecular clouds. 
Newly formed stars are associated with these dense molec- 
ular clouds and strongly influence the physical state of the 
rest of the gas in different forms through radiative and me- 
chanical inputs. The physical conditions in the multiphase 
ISM depend on the UV background radiation field, metallic- 
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ities, dust content and the den sity of cosmic rays (see Figs. 5, 
6 and 7 in IWolfire et ai1ll995l ). In addition, the filling factor 
of th e different phases depends sensitiv ely on the supernova 
rate (|de Avillez fc Breitschwerdl 120041 ). Therefore, detect- 
ing and studying the multiphase ISM in external galaxies 
has great importance for our understanding of galaxy evo- 
lution. 

Damped Lyman-a systems (DLAs) are the highest 
H I column density absorbers seen in QSO spectra, with 
N(R i)> 2 x 10 20 cm -2 . These abso rbers trace the bulk of 
the neutral hydrogen at 2 < 2 < 3 jProchaska et al.ll2005l ; 
iNoterdaeme et al]|2009h and have long been identified as re- 
vealing the interstellar medium of the high-red shift precur- 
sors of present day galaxies (for a review see, IWolfe et all 
l2005f ). 

The typical dust-to-gas ratio of DLAs, is less than 
one tenth of that observed in the local ISM, and only a 
small fraction (< 10%) of DLAs show detec t able amounts 
of molecular hydrogen JPetitiean et al.l |2000| ; iLedoux et alj 
120031 ; INoterdaeme et aljboQgt i with the dete ction rate being 
corre lated to the dust content of the gas (|Petitiean et al.l 
120061 ). The estimated temperature and molecular fraction in 
these systems are consisten t with them originating from the 
CNM (|Srianand et alj2005l ). It has been shown recently that 
strong C 1 absorbers detected in low-resolution Sloan Digital 
Sky Survey (SDSS) spectra are good candidates for H2 bear- 
ing systems. Indeed these absorbers have yielded the first 
detections of CO molecules in high-z DLA s (|Srianand et al.l 
120081 ; INoterdaeme et all [20091 . l2010l . l201ll ). The properties 
of these absorbers are similar to those of translucent molec- 
ular clouds. The fact that no DLA is found to be associated 
with a dense molecular cloud, a fundamental ingredient of 
star-formation, is most certainly related to the large extinc- 
tion that these clouds are expected to produce and/ or the 
small size of such regions (|Zwaan fc Prochaskall200q ) mak- 
ing detections difficult. 

Thus, most DLAs detected in opt ical spectroscopic sur - 
veys seem to probe the diffuse H I gas jPetitiean et al.ll2000l ). 
However, a bout 50% of the DLAs show detectable C 11* 
absorption (| Wolfe et al.ll2008l ). and IWolfe et all (|2003l ) ar- 
gued that a considerable fraction of the C 11 * absorption in 
DLA s originates from CNM gas (see however lSrianand et al.l 
2005). Detection of 21-cm absorption is the best way to 
estimate the CNM fraction of DLAs a s it is sensitive to 
both iV(H 1) and thermal state of the gas (|Kulkarni fc Heilesl 
1 19881 ). This is why it is important to search for 21-cm 
absorption in DLAs over a wide redshift range. While a 
good fraction of DLAs/sub-DLAs preselected through Mg n 

absorption seems to show 21 - cm absorption at z ~ 1 . 3 

(see for example I Gupta et al.l 12009] ; iKanekar et alj 120091 ) , 
searches for 21-cm absorption in D LAs at z a b s > 2 have 
most l y resulted in null dete ctions fsee lKanekar fc Chengalurl 
120031 ; ICurran et ail l201Ch with only four detections re- 
port^ ed_tfil_jic^seeJmifeetaD 1 19851 ; IKanekar et"aH l200fj . 
120071 ; I York et al.ll2007l ). Thelow detection rate of 21-cm 
absorption in high-z DLAs can be related to either the 
gas being warm (i.e high spin t emper ature, Ts, as sug- 
gested by IKanekar fe Chengalurl (|2003l )) and/or the low 
value of covering factor (f c ) through high-z geometric ef- 
fects (|Curran fc Webbll2006D . 

The best way to address the covering factor issue is to 
perform milliarcsecond scale spectroscopy in the redshifted 



Table 1. Log of GBT and GMRT observations to search for 21- 
cm absorption 



Source 


Tele- 


Date 


Time 


BW 


Ch. 


name 


scope 














yy-mm-dd 


(hr) 


(MHz) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


J0407-4410 


GBT 


2006-10-20 


4.7 


0.625 


512 


(CTS247) 




2007-01-05 












2007-01-06 












2007-01-08 








J0733+2721 


GBT 


2007-12-05 


10.7 


1.25 


512 






2007-12-06 








J0801+4725 


GMRT 


2006-12-22 


10.8 


1 


128 






2006-11-23 








J0852+2431 


GBT 


2009-08-07 


5.6 


1.25 


1024 






2009-08-08 












2009-09-06 








J1017+6116 


GBT 


2008-10-15 


1.5 


1.25 


1024 






2008-10-16 












2008-10-18 








J1242+3720 


GMRT 


2007-06-08 


6.1 


0.5 


128 






2007-06-10 








J1337+3152 


GMRT 


2009-01-13 


6.2 


1 


128 






2009-03-17 


7.8 


0.25 


128 


J1406+3433 


GBT 


2009-03-05 


8.0 


1.25 


1024 






2009-03-06 












2009-03-07 












2009-03-08 












2009-04-07 








J1435+5435 


GMRT 


2007-06-10 


6.7 


1 


128 



Column 1: Source name. Column 2: Telescope used for 21-cm 
absorption search. Column 3: Date of observations. Column 4: 
Total time on source i.e. after excluding telescope set-up time 
and calibration overheads. Columns 5 and 6: Spectral setup for 
the observations i.e. bandwidth (BW) and number of spectral 
channels respectively. 



21-cm line using very long baseline inte rferometry (VLBI ) 
to measure the extent of absorbing gas (|Lane et alJ|200Ch . 
Unfortunately due to limited frequency coverage and sen- 
sitivity of the receivers available with VLBI such studies 
cannot be extended to high redshift DLAs. Alternatively, 
the core fraction measured in the milliarcsecond scale im- 
ages caii_l2ej]sed to_get an estimate of the covering factor 
fsee lBriggs et alJll98a IKanekar et al.ll2009f ). Here one as- 
sumes that the absorbing gas completely covers at least the 
emission from the milliarcsec scale core. Therefore, to ad- 
dress this issue, one needs, not only to increase the number 
of systems searched for 21-cm absorption but also to per- 
form milliarcsecond scale imaging of the background radio 
sources. 

We report here the results of a search for 21-cm absorp- 
tion in 10 DLAs at z > 2 we have carried out using GBT 
and GMRT, complemented by L-band VLBA images of the 
background QSOs. This survey has resulted in the detec- 
tion of 21-cm absorption in the z a bs= 3.1745 DLA towards 
J1337+3152. A detailed analysis of this syste m and two sub- 
DLAs close to this system are presented in ISrianand et al.l 
l|2010[ ). Section [2] presents the details of our sample. In Sec- 
tion 3 we present the details of GBT and GMRT spec- 
troscopic observations, VLBA continuum observations, and 
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Figure 1. SDSS spectra showing the Lyct lines for 12 DLAs in our sample. The best fitted Voigt profiles (solid curves) together with the 
associated lcr errors (shaded regions) are over-plotted. The dotted curve gives the best fitted continuum (in some cases the continuum fit 
incl udes the emission line s also). We have used VLT UVES spectra to get 7V(H i) in the case of z a bs = 3.1745 system to wards J1337+3152 
Csee lSrianand et aI.|[2oTol) and z abs = 2.595 and 2.622 systems towards J0407-4410 (CTS 247, see lLedoux et al.ll2005^ . 
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data reduction. The detection rate of 21-cm absorption in 
DLAs is discussed in Section [4] In Section [5] we study the 
correlations between the parameters derived from 21-cm ob- 
servations, JV(H i), metallicity and redshift. In Section[6]we 
study the relation between 21-cm and H2 absorption. The 
results are summarized in Section [7] In this work we assume 
a flat Universe with Ho = 71 km s _1 Mpc -1 , fi m = 0.27 and 
Q A = 0.73. 



2 THE SAMPLE OF DLAS 

To construct our sample, we cross-correlated the overall 
sample of DLA-bearing QSO sightlines from SDSS-DR7 
(INoterdaeme et al 1120091 . including systems that are not part 
of the published statistical sample used to measure Qui) 
with the VLA FIRST catalog to identify DLAs in front of 
radio sources brighter than 50 mjy at 1.4 GHz. We excluded 
radio sources with the DLA 21-cm absorption frequencies 
redshifted into GBT and GMRT frequency ranges known 
to be affected by strong radio frequency interference (RFI) . 
There are 13 DLAs from the SDSS-DLA catalog that sat- 
isfy these conditions. In addition, there are 4 DLAs along 
the sight line towards J0407-4410 (also known as CTS 247) 
at z abs = 1.913, 2.550, 2.595 & 2.622. Two of these, at 
z a b s = 2.595 and 2.622, have redshifted 21-cm absorption 
frequency in the relatively RFI free frequency range of GBT. 
Including these two DLAs towards CTS 247, we have a sam- 
ple of 15 DLAs for which a search for 21-cm absorption 
was carried-out using either GMRT or GBT. We observed 
14 DLAs, (the exception is the Zabs— 3.079 system towards 
J1413+4505), but obtained useful spectra for only 10 DLAs. 
In addition, we have obtained milliarcsecond scale images at 
1.4 GHz for all the QSOs except CTS 247 to understand the 
role of radio structure in detectability of 21-cm absorption 
in DLAs. The details of the GBT, GMRT and VLBA obser- 
vations are given below. 

The Lyman-a profiles for 12 DLAs selected from the 
SDSS-DLA catalog are shown in the Fig. [1] The H I column 
density for each of these DLAs has been estimated using 
Voigt profile fits to the Lyman-a absorption line. The QSO 
continuum was approximated by a lower order spline using 
absorption free regions on both sides of the H I trough (dot- 
ted curves in each panel). In addition, special care was taken 
to fit the emission line profiles whenever the Lya absorption 
is close to QSO emission lines. For the remaining three DLAs 
in our sample, Za.bs=3.1745 DLA towards J1337+3152 and 
the two DL As towards CTS 247, w e use the column densitie s 
measured bv lSrianand et all (|2010T ) and lLedoux et all ||2003l ) 
respectively from high resolution VLT UVES spectra. 



3 DETAILS OF OBSERVATIONS AND DATA 
REDUCTION 

3.1 The GBT and GMRT observations 

We observed our sample of 14 DLAs using the GBT prime 
focus receivers PFl-340MHz and PFl-450MHz, and the 
GMRT P-band receiver. Although we selected DLAs such 
that the redshifted 21-cm absorption frequencies were not 
affected by strong RFI, no useful data could be obtained for 

4 absorption systems either due to RFI or other technical 



reasons. The observing log for the remaining 10 DLAs and 
the spectral set-up used for these observations are provided 
in Table [T] GBT observations were performed in the stan- 
dard position-switching mode with typically 5 min spent on- 
source and 5 min spent off-source. The data were acquired in 
the orthogonal polarization channels XX and YY. We used 
the GBT spectral processor as the backend for these observa- 
tions. The two DLAs towards CTS 247 were observed simul- 
taneously using two bands of 0.625 MHz split into 512 chan- 
nels. For the GMRT observations, typically a bandwidth of 
0.5 or 1 MHz split into 128 frequency channels was used. 
The data were acquired in the two orthogonal polarization 
channels RR and LL. For the flux density /bandpass calibra- 
tion of GMRT data, standard flux density calibrators were 
observed for 10-15 min every two hours. A phase calibrator 
was also observed for 10 min every ~45 min to get reliable 
phase solutions. 

We used NRAO's GBTIDL package to develop a 
pipeline to automatically analyse the GBT spectral-line 
data sets. After excluding time ranges for which no useful 
data were obtained, the data were processed through this 
pipeline. The pipeline calibrates each data record individ- 
ually and flags the spectral channels with deviations larger 
than 5o" as affected by RFI. After subtracting a second or- 
der baseline these data are averaged to produce baseline (i.e. 
continuum) subtracted spectra for XX and YY. The base- 
line fit and statistics for the flagging are determined using 
the spectral region that excludes the central 25% and last 
10% channels at both ends of the spectrum. If necessary, a 
first-order cubic spline was fitted to the averaged XX and 
YY spectra obtained from the pipeline, which were then 
combined to produce the Stokes-I spectrum. The spectrum 
was then shifted to the heliocentric frame. The multi-epoch 
spectra for a source were then resampled onto the same fre- 
quency scale and combined to produce the final spectrum. 

The GMRT data were reduced using the NRAO AIPS 
package foll o wing the standard procedures described in 
iGupta et all (|2006r i. Special care was taken to exclude the 
baselines and time stamps affected by RFI. The spectra at 
the quasar positions were extracted from the RR and LL 
spectral cubes and compared for consistency. If necessary, 
a first-order cubic-spline was fitted to remove the residual 
continuum from the spectra. The two polarization channels 
were then combined to get the stokes I spectrum which was 
then shifted to the heliocentric frame. 

The FWHM of the GBT beam at 400 MHz is 30' and 
the rms confusion is 500 mjy. This is comparable to the flux 
densities of the background radio sources observed with the 
GBT. Therefore, to correct for the effect of other confusing 
sources in the GBT beam and determine the QSO flux den- 
sities at the redshifted 21-cm frequency, we observed these 
with the GMRT at 610 and 325 MHz. For these observations 
we have used 32 MHz bandwidth. Details of these GMRT 
observations and the measured flux densities are provided 
in Table [2] For J1406+3433, the 325 MHz flux density is 
taken from the Westerbork Northern Sky Survey (WENSS). 
We interpolate these flux density measurements to deter- 
mine the flux densities at redshifted 21-cm frequencies for 
the quasars observed with the GBT. Since, flux densities 
for these 5 QSOs are not measured at the same epoch as 
the GBT spectroscopic observations, in principle, radio flux 
density variability can affect our estimates of 21-cm optical 
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Table 2. GMRT low- frequency flux density measurements for the DLAs observed with the GBT 



Source name 


■S610MHZ 


Date 


•S , 325MHz 


Date 




(mJy) 


yy-mm-dd 


(mJy) 


yy-mm-dd 


(1) 


(2) 


(3) 


(4) 


(5) 


J0407-4410 


124 


2006-05-23 


52 


2007-11-24 


J0733+2721 


248 


2007-11-04 


549 




J0852+2431 


198 


2008-12-26 


237 


2009-03-17 


J1017+6116 


274 




266 




J1406+3433 


165 




185+ 


(WENSS) 



Column 1: Source name. Columns 2 and 4: GMRT flux density measurements at 610 and 325 MHz respectively. Columns 3 and 5: 
Dates of the 610 and 325 MHz observations respectively, 
t From the WENSS catalog. 



Table 3. Results from the GBT and GMRT observations 



Source name 


-em 


-.1 1 >S 


logAT(H i) 


Si. 4 GHz 


5 


Spectral rms 






h 








(cm- 2 ) 


(mJy) 


(kins- 1 ) 


(mjyb- 1 eh" 1 ) 


(mJy) 


(kms- 1 ) 


(K) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


J040718-441013 


3.020 


2.595 


21.05±0.10 




3.7 


5.9 


67 


<1.61 


>382 


J040718-441013 




2.622 


20.45±0.10 






7.1 


67 


<1.93 


>81 


J073320.49+272103.5 


2.938 


2.7263 


20.25±0.20 


240 


3.8 


3.4 


451 


<0.14 


>692 


J080137.68+472528.2 


3.276 


3.2235 


20.80±0.15 


78 


7.0 


1.5 


164 


<0.22 


>1563 


J085257.12+243103.2 


3.617 


2.7902 


20.70±0.20 


160 


3.9 


3.9 


228 


<0.32 


>850 


J101725.89+611627.5 


2.805 


2.7681 


20.60±0.15 


477 


3.9 


1.2 


268 


<0.29 


>758 


J124209.81+372005.7 


3.839 


3.4135 


20.50±0.30 


662 


3.6 


3.6 


615 


<0.11 


>1567 


J133724.69+315254.5 


3.174 


3.1745 


21.36±0.10 


83 


6.9 


1.3 


69 


2.08±0.17 


600 +220 


J140653.84+343337.4 


2.566 


2.4989 


20.20±0.20 


167 


3.6 


3.0 


178 


<0.31 


>356 


J143533.78+543559.4 


3.811 


3.3032 


20.30±0.20 


96 


7.1 


1.5 


145 


<0.26 


>418 



Column 1: Source name. Column 2: QSO emission redshift. Column 3: Absorption redshift of DLAs as determined from the metal 
absorption lines. Column 4: H I column density. Column 5: Flux density at 1.4 GHz from the FIRST catalog. Columns 6 and 7: 
Spectral resolution and rms for the survey spectrum. Column 8: Continuum flux density of source at the redshifted 21-cm absorption 
frequency. Column 9: Integrated 21-cm optical depth or 3-sigma upper limit to J rdv for the equivalent spectral resolution of 
lOkms- 1 . Column 10: Ratio of spin temperature and covering factor of absorbing gas. 

depths for the corresponding DLAs. However, this effect is 
much smaller than the error caused by confusion from other 
sources in the beam and should not affect the statistical 
results derived later in the paper. 

The GBT and GMRT observations of our DLA sam- 
ple have resulted in useful 21-cm absorption spectra for 10 
DLAs. These spectra are presented in Fig. [21 GBT spectra, 
typically acquired at a spectral resolution, of ~1 km s^have 
been smoothed to ~4km s for presentation. The 21-cm 
absorption is detected only for one DLA (i.e. z a b s = 3.1745 
DLA towards J1337+3 152) and a det a iled a nalysis of this 
system is presented in ISrianand et al.l (|2010T ) . None of the 
other "absorption-like features", marked as shaded regions, 
are reproduced in spectra from different polarizations and 
epochs, but are due to RFI. For CTS247b (i.e for z abs = 
2.622 DLA towards CTS247) these features are present only 
in one polarization at certain times. For J0852+2431 and 
J1017+6116, using a combination of high spectral resolu- 
tion (~1 km s -1 ) and/or multi-epoch observations we rule 
out the possibility of these features being real 21-cm absorp- 
tion. Details of the optical depth measurements and other 
observational results for all the 10 DLAs are summarized in 
Tabled 



3.2 Continuum observations with VLBA 

The sample of quasars presented here was observed as part 
of a larger VLBA survey to obtain milliarcsecond scale im- 
ages for QSOs with foreground DLAs and Mg II systems, 
and understand the relationship between radio structure and 
detectability of 21-cm absorption. We have observed using 
VLBA 21-cm receiver for 11 hrs and 18 hrs on 21/02/2010 
and 10/06/2010 respectively. We used eight 8 MHz baseband 
channels, i.e. the total bandwidth of 64 MHz. Each baseband 
channel was split into 32 spectral points. Both the right and 
left-hand circular polarization channels were recorded. Two 
bit sampling and a post-correlation time resolution of 2 sec- 
onds were used. 

The observations were done using nodding-style phase- 
referencing with a cycle time of ~5min, i.e. 3min on the 
source and ~1.5min on the phase-referencing calibrator. 
The phase-referencing calibrators were selected from the 
VLBA calibrator survey (VCS) at 2.3 and 8.6 GHz (Table©. 
In order to improve the uv-coverage, the total observing time 
was split into snapshots over a number of different hour an- 
gles. Each source, except CTS247 which was excluded due to 
observational constraints, was typically observed for a total 
of ~30min. During both observing runs, strong fringe find- 
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CTS247a CTS247b 




395 396.8 395.4 393 393.8 393.4 



J0733+2721 J0801+4725 




405.6 406 406.4 330 330.4 

Heliocentric Frequency (MHz) 

Figure 2. GBT and GMRT spectra of DLAs in our sample. Shaded regions mark features that are due to RFI. The arrows in the case 
of non-detections indicate the expected positions of the 21-cm absorption. In the case of 2 a bs= 3.1745 system towards J1337+3152 we 
show the high resolution 21-cm absorption spectrum only. Two arrows in the case of CTS247a indicate the expected position of 21-cm 
absorption from the two H2 components. 



ers/bandpass calibrators such as J0555+3948, J0927+3902, 
J1800+3848 and J2253+1608 were also observed every ~3 hr 
for 4-5 min. 

Data were calibrated and imaged using AIPS and 
DIFMAP in a standard way. Global fringe fitting was per- 
formed on the phase-referencing calibrators. The delays, 
rates and phases estimated from these were transferred to 
the sources which were then self-calibrated until the final 
images were obtained. Radio sources were characterised by 
fitting Gaussian models to the self-calibrated visibilities. 
VLBA maps of the 13 QSOs are shown in Fig. [3] and the 
results of model fitting are listed in columns #6-12 of Ta- 
ble [321 

Non-detection of 21-cm absorption in a DLA could be 



due to the small covering factor of the absorbing gas. The 
typical spatial resolution achieved in our VLBA observations 
is ~8 mas. If the extent of absorbing gas is of the order of the 
scales probed by our VLBA observations (i.e. >20pc) then 
we expect the detectability of 21-cm absorption to depend 
on the fraction and spatial extent of radio flux density de- 
tected in these images. In column # 14 of Table ET2"1 we give 
the ratio of total flux densities detected in the VLBA and 
FIRST images at 20cm, i.e. /vlba- The last column of this 
table gives the largest linear size (LLS), i.e the separation 
between the farthest radio components, of the radio source 
at the redshift of the DLA. Out of the 13 QSOs presented 
in Fig. [3 that have DLAs along their line of sight, we have 
21-cm absorption spectra for only 8 DLAs. For the DLA to- 
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Figure 3. Contour plots of VLBA images at 1.4 GHz. The rms in the images are listed in Table [3721 At the bottom of each image the 
restoring beam is shown as an ellipse, and the first contour level (CL) in mjybeam - 1 and FWHM are noted. The contour levels are 
plotted as CLx(— 1, 1, 2, 4, 8,...) mjybeam" 1 . Depending upon the detailed structure of the radio sources, the emission could be more 
extended at the redshifted 21-cm frequencies. 



Table 5: Results from the VLBA data. 



Source 


■^abs 


Right 


Declination 


rms 


Comp. 


S 


r 


6 


a 


b/a 




S T 


/VLBA 


LLS 


name 




ascension 






























(J2000) 


f ~X c\r\r\r\\ 

(J 2000) 


/ T I — 1 \ 

(mjy beam J 




(mJy) 


(mas) 


( ) 


(mas) 




( ) 


(mJy) 




(PC) 


1 1 \ 
(!) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


J0733+2721 


2.7263 


07 33 20.4830 


+27 21 03.430 


0.3 


1 


97 





■ 


4.52 


0.05 


-88 


240 


0.62 


348 












2 


3 


9.1 


-129 


13.24 


0.00 


-18 


















3 


19 


32.7 


-140 


5.98 


0.26 


32 


















4 


29 


10.8 


51 


8.37 


0.44 


-17 








J0801+4725 


3.2235 


08 01 37.6930 


+47 25 28.082 


0.2 


1 


28 





- 


6.77 


0.00 


72 


78 


0.67 


342 












2 


8 


9.34 


-112 


7.09 


0.00 


72 


















3 


16 


44.7 


-112 


8.63 


0.22 


64 








J0816+4823 




08 16 19.0044 


+48 23 28.490 


0.2 


1 


42 





- 


4.59 


0.24 


87 


69 


0.75 


68 












2 


10 


9.0 


132 


8.75 


0.43 


-63 








J0839+2002 




08 39 10.8970 


+20 02 07.391 


0.2 


1 


113 





- 


4.91 


0.24 


84 


130 


0.87 


<39 


J0852+2431 


2.7902 


08 52 57.1211 


+24 31 03.271 


0.2 


1 


78 





- 


4.92 


0.31 


71 


160 


0.55 


175 












2 


10 


21.9 


42 


15.2 


0.16 


53 








J1017+6116 


2.7681 


10 17 25.8865 


+61 16 27.414 


0.5 


1 


388 





" 


1.50 


0.68 


55 


477 


0.86 


38 












2 


24 


4.7 


145 


2.45 


0.00 


70 








J1223+5037 




12 23 43.1740 


+50 37 53.344 


0.5 


1 


96 





- 


2.50 


0.00 


78 


229 


0.60 


554 












2 


16 


13.6 


80 


3.01 


0.83 


-19 


















.3 


25 


71.4 


78 


8.61 


0.00 


89 








J1237+4708 




12 37 17.4413 


+47 08 06.964 


0.2 


1 


64 







3.17 


0.20 


-74 


80 


0.80 


< 27 


T1 0/1 J_^79n 


0.4100 


10/19 HQ £1 01 

±Z 4Z uy.oizi 


j_q*7 on ^qo 
-\-o( zu uo.oyz 


U-U 


1 

1 


o4o 


n 
U 






U. ( D 




zz 


DDZ 


1 nn 


\ 14 


J1337+3152 


3.1747 


13 37 24.6931 


+31 52 54.642 


0.2 


1 


83 







3.85 


0.38 


74 


83 


1.00 


< 30 


J1406+3433 


2.4989 


14 06 53.8532 


+34 33 37.339 


0.4 


1 


127 







3.24 


0.22 


-23 


167 


0.87 


153 












2 


18 


18.7 


-30 


23.79 


0.23 


-25 








J1413+4505 




14 13 18.8652 


+45 05 22.990 


0.2 


1 


105 







2.19 


0.42 


-67 


140 


0.88 


216 












2 


12 


3.4 


-77 


5.37 


0.00 


-86 


















3 


6 


27.9 


-72 


6.28 


0.03 


-67 








J1435+5435 


3.3032 


14 35 33.7812 


+54 35 59.312 


0.2 


1 


31 







2.66 


0.17 


-29 


96 


0.55 


155 












2 


17 


20.4 


155 


4.00 


0.47 


-36 


















3 


5 


7.6 


153 


5.52 


0.00 


-7 









<u Column 1: Source name. Column 2: absorption redshift. Columns 3 and 4: right ascension and declination of component-1 (see column 6) from the multiple Gaussian fit 

H to the source, respectively. Column 5: rms in the map in mjybeam -1 . Column 6: component id. Column 7: flux density of the component in mJy. Columns 8 and 

C3 9: radius and position angle of the component with respect to component-1, respectively. Columns 10, 11 and 12: major axis, axial ratio and position angle of the 

§ deconvolved Gaussian component, respectively. Column 13: flux density in mJy from FIRST/NVSS. Column 14: c/ is the ratio of 1.4 GHz flux density in VLBA image 

£ to that in the FIRST image, Column 15: largest projected linear size in pc. 



CO 



21-cm absorption in DLAs 9 



Table 4. Details of phase-referencing calibrators used for the 
VLBA observations 



100 



Source 


Calibrator 


Separation 






(degrees) 


(1) 


(2) 


(3) 


J0733+2721 


J0732+2548 


1.5 


J0801+4725 


J0754+4823 


1.5 


J0816+4823 


J0808+4950 


1.9 


J0839+2002 


J0842+1835 


1.6 


J0852+2431 


J0856+2111 


3.5 


J1017+6116 


J1031+6020 


2.0 


J1223+5037 


J1227+4932 


1.3 


J1237+4708 


J 1234+4753 


0.9 


J1242+3720 


J1242+3751 


0.5 


J1337+3152 


J1329+3154 


1.6 


J1406+3433 


J 1416+3444 


1.9 


J1413+4505 


J1417+4607 


1.2 


J1435+5435 


J1429+5406 


1.0 



Column 1: Source name. Column 2: Phase-referencing 
calibrator. Column 3: Separation between the radio source and 
phase-referencing calibrator. 



J1435+5435 




, " CL = O.f 

1:0.013" xo.oo5" : -is 
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Figure 3. Continued. 

wards J08 16+4823, we use th e 21-cm absorption measure- 
ment from lCurran et all (|201(t ). Thus we have a sample of 9 
DLAs with both 21-cm absorption measurements and VLBA 
21-cm maps for the background QSOs. The /vlba for this 
sample ranges from 0.6 to 1, and LLS from <15 pc to 340 pc. 

In the absence of VLBI spectroscopy at the redshifted 
21-cm line frequency, the ratio of VLBA core flux density to 
the total flux density measured in the arcsecond scale im- 
ages (called core fraction c A has been used as an indicator of 
the c overing factor f c fsee lBriggs et al] fl989: Ka nekar et al.l 
I2009T ). Here we use the term 'core' to refer to the flat spec- 
trum unresolved radio component coincident with the opti- 
cal QSO in the VLBA image. 

For J1242+3720 and J1337+3152 the radio source is 
modelled as a single unresolved component and within the 
uncertainties all the flux density in the arcsecond scale 
FIRST images is recovered in our VLBA images (Fig. [3]). 
Both of these sources, have flat spectra suggesting the radio 
emission even at lower frequencies originates predominantly 
from the compact core. Therefore we take / c =l for this case. 

The radio source J1017+6116 has an inverted radio 




1000 



2000 



3000 



T[ (K) 



Figure 4. The allowed range of fraction of DLAs having har- 
monic mean spin temperature Ts greater than a limiting value 
T| as a function of Tg. We use only those DLAs for which Cf 
measurements are available. The lower envelope of the shaded 
region is obtained considering all the lower limits on Ts as mea- 
surements. The upper envelope is obtained assuming all the lower 
limits as measurements with Tg > Tl. 



spectrum (a=— 0.4) and 86% of the flux density in the 
FIRST image is recovered in the VLBA image, 94% of which 
is contained in the main unresolved component (Fig. [3}. For 
another 3 flat-spectrum sources with 21-cm absorption mea- 
surements, i.e. J0816+4823, J0852+2431, and J1406+3433, 
more than 80% of the VLBA flux density is present in a sin- 
gle unresolved component. For these 4 sources, based on the 
flat spectral index, the dominant component in the VLBA 
image can be identified with radio core/optical QSO. There- 
fore for the 6 sources mentioned above, we have estimated 
the core fraction, Cf, and used it as the covering factor, 
f c , of the gas. The remaining three sources, J0733+2721, 
J0801+4725 and J1435+5435, exhibit multiple components 
in their VLBA 20-cm images. The identification of the com- 
ponent coincident with the QSO, and the estimation of their 
Cf from VLBA images for these three sources are highly un- 
certain. 



4 DETECTABILITY OF 21-CM ABSORPTION 

In this Section we investigate the detectability of 21-cm ab- 
sorption in DLAs and the implication of non-detections for 
the physical state of the H I gas. It is clear from the last 
column of Table [3] that for most of the DLAs, our data has 
good sensitivity to detect Ts// C ~100 K gas. 

The H I 21-cm absorption is detected only in the 
z abs = 3.1745 system towards SDSS J1337+3152. This is 
one of the weakest radio sources in our sample (with a 3a 
J rdv limit of 0.4 km s _1 ). However, thanks to high N(H. i) 
our spectrum is sensitive enough to detect any gas with 
Ts < 3100 K. This source is unresolved in our VLBA ob- 
servations (see Fig [3}- The L-band flux density measured in 
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our VLBA image is consistent with that measured by the 
FIRST survey. Therefore, the core fraction is, a ~ 1, and 
the size of the VLBA beam is less than 30 pc at the redshift 
of the absorber. The spin-temperature, measured from the 
ratio of 21-cm optical depth and the iV(H i) column den- 
sity derived from the Lyman-a trough, is 600ljgg K which 
is consistent with the upper limit on Ts obtained from the 
width of th e single component Gaussian fit to the 21-cm 
absorption (|Srianand et alJboiCh . 

In Table [S] we provide various details of our measure- 
ments together with the previous measurements at z > 2 
from the literature. We present the results dividing the sam- 
ple into three groups. These are systems with 21-cm detec- 
tions (five systems), systems with 21-cm absorption upper 
limits with (twelve systems) and without (eleven systems) 
high-resolution optical spectra from which to derive accu- 
rate metallicities. The first two groups are used to investi- 
gate the connection between UV measurements and 21-cm 
optical depth. In all cases the 3<r upper limits on the inte- 
grated 21-cm optical depth are computed assuming a line 
width of 10 km s _1 . 

The 21-cm detection rate from our sample, without 
putting any sensitivity limit, is 10%. This is 13% when we 
restrict to J rdv limit of 0.4 km s _1 ( the limit achieved in 
the case of J1337+3152 where we have 21-cm detection). 
Taken at face value, the extended sample listed in Table [6] 
gives a 21% detection rate for J rdv limit of 0.4 km s~ . 
For a J rdv limit of 0.2 km s~ we get the detection rate of 
28%. However, these may not be representative values as the 
list of systems compiled from the literature may be biased 
towards detections as some authors may not have reported 
their non-detections systematically 

Since we know N(H i) from the damped Lyman-a line, 
the detection limit on the integrated optical depth implies a 
lower limit on the ratio Ts// C . The Tg// C measurements are 
reported in column 7 of Table[6] In column 6 of this table, we 
give the core fraction c/. As mentioned above, c/ is basically 
the ratio of flux density in the unresolved core seen in VLBA 
images to the total flux density measured in the arcsecond 
scale FIRST images. For the objects from the literature we 
use the c/ values given in lKanekar et al] |2009l ). These mea- 
surements were made at 327 MHz, close to the re dshifted 
21-cm frequencies. Following iKanekar et al.l l|2009l ) we use 
core fraction (c/) as the estimate of the covering factor (/ c ). 
The Ts measurements given in column 8 of Table are ob- 
tained by assuming f c — Cf. 

In Fig. [4] we plot the percentage of DLAs having Tg 
greater than a limiting value T s as a function of Tg for sys- 
tems with Ts measurements given in column 8 of Table [6] 
The lower envelope of the shaded region is obtained con- 
sidering all the lower limits on Ts as measurements. The 
upper envelope is obtained assuming all the lower limits as 
measurements with Ts > Tg. It is clear from the figure that 
more than 50% of the DLAs have T s > 700J\ . Remember 
that the Tg measured in an individual DLA is the harmonic 
mean temperature of different phases that contribute to the 
observed iV(H i). Assuming that the gas is simply a two 
phase medium with similar covering factors the fraction of 
H I in the CNM (called /(CNM)) can be written as, 




/(CNM) 



rpW 



T K 



T s c 



(1) 



Figure 5. Metallicity vs Ts/fc- The vertical dashed line marks 
the median metallicity measured in our sample. 



where, T s c and T^ are the spin-temperature of the CNM 
and WNM respectively. ISrianand et al.l (|2005l ) have noticed 
that the H I phase traced by the H2 absorption has tem- 
perature typically in the range 100-200 K. Thus we consider 
the CNM temperature to be 200 K (instead of 70 K as seen 
in CNM of the Galaxy) so that the /(CNM) we get will 
be a conservative upper limit. Assuming T$ ~ 200 K and 
T^ ~ 10 4 K, T s = 700 K can be obtained for a combi- 
nation of /(CNM) = 0.27 and /(WNM) = 0.73. Therefore 
/(CNM) is less than 0.27 in at least 50% of the DLAs. Note 
that choosing Tg* ~ 8000 K (as suggested for the Galactic 
ISM) instead of the 10 4 K used here, does not change the 
results appreciably. 

We estimate /(CNM) for the four 21-cm detections (ex- 
cluding J0501-0159 (B0458-020) for which we do not have 
the covering factor value). Apar t from J0314+43 14 (3C082) 
which seems to be a special case l|York et al.ll2007h . the CNM 
seems to represent roughly 20 to 30% of the total iV(H 1) 
measured in these DLAs. For individual non-detections, we 
can calculate conservative upper limits of the fraction of 
A(H 1) in the CNM phase assuming T s c = 200 K. The val- 
ues of /(CNM) are given in column #9 of Table [6] for sys- 
tems with f c measurements. The upper limits vary between 
0.10 and 1.0 with a median value of 0.23. Thus the analysis 
presented here, under the assumption that f c = Cf, sug- 
gests that most of the neutral hydrogen in high-z DLAs is 
warm. This is very m uch consistent with the conclusion of 
iPetitiean et~aT] ||2000h based on the lack of H2 detections in 
most high-z DLAs. 



5 RESULTS OF CORRELATION ANALYSIS 

In this Section we explore correlations between the 21-cm 
optical depth and other observable parameters. As we have 
only a few 21-cm detections and mostly upper limits we 
use survival analysis and in particular the generalised rank 
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Table 6. Summary of 21-cm searches in z > 2 DLAs. Column 1: QSO. Column 2: absorption redshitft. Column 3: log iV(H i). Column 
4: Integrated optical depth. Column 5: Reference for J rdv given in column 4. Column 6: the core fraction Cf. Column 9: the fraction of 
CNM (see the text for its definition). Column 10: the H2 fraction and Column 11: References for JV(H 1) and/or f(H2) measurements. 



QSO 




logA r (H 1) 


J rdv 


Refs 


c f 


Ts/fc 


T S 


/(CNM) 


log f(H 2 ) 


Refs 








(km s _1 ) 






(K) 


(K) 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 








21-cm detections 












J020346.6+113445 


3.38714 


21.26+ 0.08 


0.71+0.02 


1 


0.76 


1397 


1062 


~0.19 


-6.2,-4.6 


a 


J031443. 6+431405 


2.28977 


20.30+ 0.11 


0.82+0.09 


2 




133 




~1.00 




d 


J044017.2-433309 


2.34747 


20.78+ 0.10 


0.22+0.03 


3 


0.59 


1493 


881 


~0.23 




c 


J050112.8-015914 


2.03955 


21.70+ 0.10 


7.02+0.16 


1 




390 






< -6.40 


b 


J133724.6+315254 


3.17447 


21.36+ 0.10 


2.08+0.17 


5 


1.00 


600 


600 


~0.33 


-7.00 


c 




21- 


cm non-detections having metallicity 


measurements 








J033755. 7-120412 


3.1799 


20.65+ 0.10 


<0.06 


6 


0.62 


>4057 


>2515 


< 0.08 


< -5.10 


b 


J033901. 0-013318 


3.0619 


21.10+ 0.10 


<0.06 


6 


0.68 


>11435 


>7775 


< 0.03 


< -6.90 


b 


J040733. 9-330346 


2.569 


20.60+ 0.10 


<0.12 


7 


0.44 


>1807 


>795 


< 0.25 




c 


J040718. 0-441013 


2.59475 


21.05+ 0.10 


<1.61 


8 




>380 






-2 61+ ' 17 


b 


J040718. 0-441013 


2.62140 


20.45+ 0.10 


<1.93 


8 




>80 






< -6.20 


b 


J043404.3-435550 


2.30197 


20.95+ 0.10 


<0.33 


7 




>1471 






< -5.15 


b 


J053007.9-250330 


2.81115 


21.35+ 0.07 


<0.58 


9t 


0.94 


>2103 


>1977 


< 0.10 


9 oo+0. 18 


b 


J091551. 7+000713 


2.7434 


20.74+ 0.10 


<0.37 


7 




>809 






c 


J135646. 8-110129 


2.96680 


20.80+ 0.10 


<0.33 


6 




>1042 






< -6.75 


b 


J135706. 1-174402 


2.77990 


20.30+ 0.15 


<0.14 


7 




>777 






< -5.99 


a 


J142107.7-064356 


3.44828 


20.50+ 0.10 


<0.14 


7 


0.69 


>1231 


>849 


< 0.24 


< -5.69 


b 


J234451. 2+343348 


2.90910 


21.11+ 0.10 


<0.21 


6 


0.71 


>3343 


>2373 


< 0.08 


< -6.19 


a 




21-cm non-detections having 


no metallicity measurements 








J053954.3-283956 


2.9742 


20.30+ 0.11 


<0.06 


6 


0.47 


>1812 


> 851 


< 0.23 




c 


J081618. 9+482328 


3.4358 


20.80+ 0.20 


<1.43 


9 


0.60 


>240 


> 144 


< 1.00 




a 


J073320.4+272103 


2.7263 


20.25+ 0.20 


<0.14 


8 




>692 








a 


J080137.6+472528 


3.2235 


20.80+ 0.15 


<0.22 


8 




>1563 








a 


J085257.1+243103 


2.7902 


20.70+ 0.20 


<0.32 


8 


0.49 


>854 


> 418 


< 0.48 




a 


J101725. 8+611627 


2.7681 


20.60+ 0.15 


<0.29 


8 


0.81 


>748 


> 606 


< 0.33 




a 


J124209. 8+372005 


3.4135 


20.50+ 0.30 


<0.11 


8 


1.00 


>1566 


> 1566 


< 0.13 




a 


J140501. 1+041536 


2.708 


21.07+ 0.24 


<0.19 


9 




>3369 








f 


J140501. 1+041536 


2.485 


20.20+ 0.20 


<0.08 


9 




>1080 








f 


J140653. 8+343337 


2.4989 


20.20+ 0.20 


<0.31 


8 


0.76 


>279 


> 212 


< 0.94 




a 


J143533. 7+543559 


3.3032 


20.30+ 0.20 


<0.26 


8 




>418 








a 



References in column #5: ltlKanekar et al.l ll2007fl. 2) [York et al.l J2007I'). 3) iKanekar et al.l j2006l l. 4~l lBriggs et al.l jl989T), 5 ) 
ISrianand et al ] l|2010h . 6) lKanekar fc Chengalurl l|2003h . 7) iKanekaretal ] l|2009h . 81 This paper, and 9) ICurrari^tal .1 |201Ch . t Archival 

data from GBT08A_003 (PI: Curran) was processed through our pipeline. See text f or details. 

Refere nce s for jV(H i) and/or f(H 2) measurements (colu mn # 11): a) Th i s pape r, b) rNotcrdac m"e"et al.l l|2008h . c) ISrianand et all 
j201Ch . dl lEUison et al.l 11200811 , el lAkerman et al.l ||2005|1 fl lCurran et al.l <2010h . 



correlation test ()lsobe et al.lll986i ). For this purpose we use 
the Astronomical SURVival analysis (ASURV) package. 



5.1 Metallicity vs T s // C 

Firstly we study the importance of the metallicity of the gas. 
Only 3 DLAs in our sample have measurements of metal- 
licity from high resolution optical spectroscopy. In the ex- 
tended sample at z > 2 (see Table [6} there are 17 systems 
with metallicity measurements and 21-cm spectra. In Fig.[5l 
we plot Ts/fc versus metallicity. The vertical long-dashed 
line marks the median metallicity of the points plotted in the 
figure. The only detection found in the low metallicity half 
is for z a bs = 3.1745 towards J1337+3152 reported from our 
survey. The other four detections are from the high metallic- 
ity half. The non-parametric generalized Kendall rank cor- 
relation test suggests only a weak correlation between Z 
and Ts/fc (at the 1.42<r level) with the probability that 
it can arise due to chance being 0.15. The significance is 



even lower (i.e 0.9cr with a chance probability of 0.37) when 
we use Ts (instead of Ts/f c ) for cases where we have esti- 
mated f c measurements. We wish to point out that a corre- 
lation between T a /f c and metallicity is reported in the liter- 
ature (|Curran et al ] |2007l : lKanekar et al . 2009; C urran et all 
l2010t ). The lack of correlation in our sample (with systems 
in a restricted redshift range) could either reflect redshift 
evolution of the relationship or small range in metallicity 
covered by the sample. Metallicity measurements for the re- 
maining 11 systems in Table [S] would allow us to address 
this issue in a statistically significant manner. 

We also looked at the possible correlation between the 
21-cm optical depth and the velocity width (Av) of the low 
ionization lines. This information is available for 14 sources. 
Again we find no statistically significant correlation between 
the two. This is inconsistent with the 2.2 to 2.8 a level 
weak correlat i on be tween VK(Mg n) and Av reported by 
ICurran et all ll2007l). H owever , this is consistent with the 
finding of iGupta et all (|2009h . that W(Mg n) and 21-cm 
optical depth are not correlated for a sample of 33 strong 
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Figure 6. Properties of 21-cm absorptions vs redshift and ./V(H i). The vertical dashed lines give the median value of the quantity plotted 
in the x-axis. 



Mr i i systems at 1.10< z < 1.45 (see also iKanekar et al.l 
l2009h . 



5.2 Redshift dependence 

In the left hand side panels of Fig. [6]we plot T s /f c and inte- 
grated 21-cm optical depth vs redshift. No clear correlation 
is evident in this figure. The non-parametric Kendall test 
finds no significant correlation between J rdv (or Ts) and 
z. Note our sample probes only a restricted redshift range 
in terms of time interval probed. However, the lack of cor- 
relation found here is consistent with the near constancy o f 
T s / f c as a function of redshift found bv lCurran et al.l(|2O10D . 
Understanding the redshift dependence of Ts is very impor- 
tant in particular to address whether there is any evolu- 



tion in T„ (IKanekar fc Chengalurll2003r i or geometric effects 
|Curran fc Webbll2006t ). To make an unbiased comparison 
we need to have 21-cm measurements at low z for a well 
defined sample of DLAs detected based on Lya absorption. 



5.3 Dependence on iV(H i) 

Recently ICurran et al.l |2010l ) have found a 3a level corre- 
lation between iV(H i) and T„/f c . To check whether this 
correlation holds at z > 2, we plot, in the top panels of 
Fig. [6l the integrated 21-cm optical depth as a function of 
redshift and iV(H i). We note that there is a tendency for 
more 21-cm detections in DLAs with higher JV(H i). How- 
ever, the non-parametric Kendall test finds no significant 
correlation between J rdv and iV(H i). In the bottom right 
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panel we plot Ts//c against log JV(H i). The Kendall test 
does not show any significant relation between the two quan- 
tities (1.28a with a probability of 0.2 for this to be due to 
chance). Thus we do not find any evidence for the 21-cm 
optical depth to depend on N(H i) in our sample. 



6 21-CM ABSORPTION AND H 2 

As 21-cm absorption and H2 molecules can give comple- 
mentary information on the physical state of the gas. In 
this Section, we study the relationship between these two 
indicators. There are 13 DLAs in our extended sample for 
which the expected optical wavelength range of redshifted 
H2 absorptions has been observed at high spectral reso- 
lutio n. Nine of these sources a re part of the UVES sam- 
ple of lNoterdaeme et all (|2008MSrianand et all (|2010h have 
reported the detection of H2 in J1337+3152 and here we 
report the search for H2 in the remaining three DLAs 
(z abs = 3.3871 towards J0203+1134, z abs = 2.7799 towards 
J1357-1744 and z abs = 2.9091 towards J2344+3433) . In the 
10th column of Table [5J we summarize the molecular frac- 
tion /(H 2 ) = 2N(U 2 ) / (2N(R 2 )+N(R 1))] derived for these 
13 systems. 

In 8 systems, neither 21-cm absorption nor H2 molecules 
are detected with typical upper limits of the order of 10~ 6 
for /(H2). Apart from the system at z abs — 2.6214 towards 
J 0407— 4410, the lower limits on Ts// C for the remaining 7 
systems are higher than 700 K. There are 4 cases where / c 
measurements are available. In three cases (z abs = 3.1799 to- 
wards J0337-1204, ,z abs = 3.0619 towards J0339-0133 and 
z abs = 2.9019 towards J2344+3433), the lower limit on T s 
is more than 2000 K. Th ese are in line with the suggestion 
bv lPetitiean et all (|2000l ) that the absence of H2 in most of 
the DLAs is due to the low density and high temperature of 
the gas. 

In two cases (z abs = 2.5947 towards J0407-4410 (CTS 
247) and z abs = 2.8112 towards J0530-2503 (PKS 0528- 
250)), strong H2 absorption is detected with rotational ex- 
citations consistent with the H2-bearing gas being a CNM. 
However, 21-cm absorption is not detected in either case. 
We discuss these two systems in detail below. 

Among the five 21-cm absorbers, high resolution UVES 
spectra covering the expected wavelength range of H2 
absorption are available for four systems. The excep- 
tion is the z abs = 2.28977 system towards J0314+4314 
(B0311+430). For the z abs = 2.3474 system towards 
J0440-4333 (B0438-436) the continuum flux in the ex- 
pected wavelength range is removed by high ionization lines 
from an associated system, as well as by a high-z Lyman 
limit system present along the line of sight. Below we dis- 
cuss the five systems where simultaneous analysis of H2 and 
21-cm absorption is possible. 



6.1 z abs = 3.3868 DLA towards J0203+1134 (PKS 
0201 + 113) 

Searches for 21-c m absorption in this system have yielde d 
conflicting results |de Bruvn et al.|[l99r3 ; lBriggs et al.ll 19971 ). 
Based on GMRT sp ectra taken at three different epochs, 
iKanekar et al.l (|2007t ) reported the detection of 21-cm ab- 
sorption in two components at z abs = 3.387144(17) and 



3.386141(45). Using JV(Hl)~ (1.8±0.3)xl0 21 cm" 2 they ob- 
tained T s = [955 ± 160l(.fc/0.69) K. U sing high resolu- 
tion optical spectrum. lEUison et al.l (|200ll ) have found a gas 
phase metallicity of 1/20 of solar with very little dust de- 
pletion. The gas cooling rate, log l c = —26.67 ± 0.10 erg 
s _1 Hz -1 , derived using C n* absorption is consistent with 
this DLA being p art of high-cool population defined by 
IWolfe et all |20Q3). From IKanekar et all (|2007T ) we can see 
that the strongest 21-cm absorption does not correspond to 
the strongest velocity component in either C 11* or Fe II. 

Here we report the detection of H2 absorption from J=0 
and J=l levels originating from both Lyman and Werner 
bands (see Fig. [7|). A single component Voigt profile fit re- 
produces the data well. As the Lyman-a forest is dense and 
the spectral signal-to-noise ratio is not very high due to 
the faintness of the QSO, we considered a range of b val- 
ues (i.e between 1 and 5 km s _1 ) to get the best fit values 
of log[iV(H 2 , J=0)] in the range 16.10-14.48 and log[AT(H 2 , 
J=l)] = 16.03—14.57. We estimated the kinetic temperature 
using the ortho-to-para ratio (i.e Tbi) and found it to be in 
the range 48—108 K for the range of b parameters considered 
above. We note that for b parameters greater than 2 km s _ , 
the H2 lines are mainly in the linear portion of the curve of 
growth and the column density estimate is insensitive to the 
assumed value. The average molecular fraction, log /(H2), 
in the range, -4.6 < log /(H 2 ) < -6.2. 

Despite the gas being cold, there is no 21-cm absorp- 
tion detected at the position of the H2 component (at z = 
3.38679) which is well separated from the 21 -cm absorption 
compo nents If we use f c = 0.76, as found by IKanekar et akl 
(;2009) using 326 MHz observations, we find log AT(H 1) 
< 19.12. This is less than 1% of the total H I column density 
measured in this system. 

Unlike most of the strong H2 systems, this system does 
not show detectable C I absorption. This means we do not 
have, unfortunately, an independent estimate of the density 
from fine-structure excitation. 



6.2 z abs = 2.5948 towards J0407-4410 (CTS 247) 

As the radio source is faint, our GBT spectrum only gives 
a weak limit on the spin temperatu re, Ts > 380 K when 
we use a line width of 10 km s" - 1 . ISrianand et "all (|2005h 
have reported log N(C n*) = 13.66±0.13. This, together 
with log JV(H 1) = 21.05±0.10, gives a gas cooling rate of 
log l c = 26.92 ±0.16. This is very close to the value lf H that 
seems to d emarcate be t ween the high and low cool systems 
defi ned by I Wolfe et al.l (120081 ) . 

iLedoux et al.l |2003l ) reported the detection of H2 
from this system. The H2 absorption is well fitted with 
two components at z abs = 2.59471 and 2.49486 with 
log N CR2) = 18.14 and 15.51 respectively (|Srianand et al.l 
120051 ). These components have T 01 = 121±10 and 91±6 K 
respectively. The average molecular fraction, log /(H2), 
is found to be — 2.421q'j2 with an average metallicity of 
— 1.0 2 ± 0.12 and moderate dust depletion |Ledoux et al.l 
I2003D . 

The absence of 21-cm absorption from this system is 
intriguing as H2 components have T~100 K. With the same 
b parameters as used to fit the H2 lines and the rms from 
the GBT spectrum, we get a 3a upper limit of J rdv = 
0.88 km s . This translates to a constraint, / c xiV(H 1) < 
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Figure 7. Voigt profile fits to H2 Lyman and Werner band absorption lines in the z^^ = 3.3868 DLA system towards J0203+1134 
(PKS 0201+113). The zero of the velocity scale is defined a t z = 3.38716. The t wo vertical lines at v = and —68 km s~ 1 show the 
locations of two 21-cm absorption components reported by llKanekar et al.ll2007h . The vertical dotted line indicates the location of H2 
absorption. In each panel we also show the error spectrum with dotted curves. 



2xl0 20 cm" 2 in the H2 components where we have assumed 
Ts — Tqi. Unfortunately we do not have a VLBA image of 
this source and it is difficult to constrain the covering factor 
of the gas. If we assume f c ~ 1 then the upper limit on 
JV(H 1) implies that the H2 component is a sub-DLA with 
log /(Ha) > -1.85. 

From the column densities of the C I fine-structure lines, 
ISrianand et al.l |2005h have constrained the particle density 
in the gas to be in the range 4.5 < hh(ciii" 3 ) < 57.3. For 
71h = 4.5 cm" 3 and f c = 1, we estimate the thickness of the 
H2 cloud (i.e A r (H i)/tih ) along the line of sight to be < 15 
pc. 

6.3 z abs = 2.0395 towards J0501 0159 (B0458-020) 

21- cm absorpti o n in t wo velocity components was reported 
by I Wolfe et all (|l985l ). The background radio source shows 
structure over a wide range of scales and t he absorbing 
cloud seems to cover most of these components (|Briggs et al.l 
1 19891 ). The estimated extent of the H I absorber is ~8 kpc 
and the spin-temperature of the system is 390 K. 

C n* absorption is det ected and the mea sured cooling 
rate is l c = -26.41 ± 0.10 \ Wolfe et alj|2008h . This is con- 
sistent with the high cool population. H2 molecules are not 



detected with an upper limit of the molecular fraction of 
log /(H2) < —6.52. This is one of the rare DLAs to show 
Lyman-a emission in the middle of the DLA absorption 
|M0ller et alJl2004lHeinmuIler et alj|2006f ). 

Based on the star format i on ra te derived from the 
Lyman-a flux. lHeinmiiller et al.l |2006l ) argued that the am- 
bient radiation field is 10 times higher than the Galactic 
UV background. This excess radiation could be the main 
reason for the absence of H2 in the gas. This picture is also 
confirmed by the lack of associated C I absorption. 



6.4 2 abs = 2.8111 towards J0530-2503 (PKS 
0528-250) 

ICarilli et ail (|l996f ) searched for 21 -cm absorption in this 
system using the WSRT and did not fin d any significant ab- 
sorption. Recently ICurran et al.l (|2010h have reported GBT 
observations of this source with a better signal-to-noise. 
They obtain a lower limit of the spin temperature of 700 K 
for f c = 1, again assuming a width of 10 km s~ . This source 
was found to be compact in the low f requency VLBA images 
with f c = 0.94 |Kanekar et al.ll2009l ). 

No measurement of cooling rate is possible for this sys- 
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tem as the the C II* absorption is blended with other satu- 
rated lines. 

This system has log iV(H i) = 21.35±0.07, 
Z = —0.91 ± 0.07 and z a b s higher than z em . This 
is the first high reds hift DLA where H2 abs o rption 
has been detected (jLevshakov fc Varshalovichl 1 19851 ; 
ISrianand fc Petitieanl 1 19981 ; ISrianand et all 120051 ). The 
high resolution UVES spectrum reveals two strong H2 
components at z abs = 2.81100 and 2.81112 with log jV(H 2 ) 
= 17.93lo.2o and 17-90tg;i4, and T01 = 167±7 and 138±12 
K respectively. 

We re-reduced the GBT data of ICurran et al.l (|2010l ) in 
the same way as our GBT data and found the rms in the 
frequency range expected for the 21-cm absorption associ- 
ated with the H2 components to be 6.7 mjy. This, together 
with the temperature Tbi and the H2 b parameter give a 
constraint of log iV(H 1) < 20 in the H2 components. 

C 1 absorption is detected only in one of the H 2 com- 
ponents (i.e at z a b s = 2.81112) (|Srianand et al.ll2005h . From 
the C 1 level populations and for the above temperature, 
we derive that the particle density is in the range 25 < 
n H (cm" 3 ) < 270. If we use the AT(H 1) limit of 10 20 cm" 2 
derived from the 21-cm absorption, and use the constraint 
we get for the hydrogen density in the case of the z a bs = 
2.81112 component, we get a limit on the thickness of the 
gas along the line of sight to be 1.33 pc. 



6.5 z abs = 3.1745 towards J1337+3152 

D etection of 21-cm and H 2 in this system was reported 
111 ISrianand et~afl (|20ld) . The weak H2 absorption is well 
aligned with the strongest metal line component but is 
slightly shifted by 2.5 km s _1 with respect to the 21-cm ab- 
sorption component. This again suggests that the two tran- 
sitions do not originate exactly from the same gas and that 
the absorptions arise from an inhomogeneous absorbing re- 
gion. 

As the H2 column density is very low and the ra- 
tio of iV(H2,J=0) to Af(H2,J=l) is close to the maxi- 
mum value allowed by the Boltzmann distribution, we 
cannot constrain the kinetic temperature of the gas. 
C n* absorption is detected with a total column density 
log N{G 11*) = 13.61±0.08. Using log AT(H 1 ) = 21.36±0.10 
we derive a gas cooling rate of log l c = —27.28 ± 0.13. Thus 
the system belongs to the low cool part of the DLA pop- 
ulation. The quasar is unresolved in our 1420 MHz VLBA 
image with a total flux density consistent with the measure- 
ment from FIRST observations. This suggests that most of 
the emission is in the unresolved component. Using Gaus- 
sian fits we estimate the limit on the largest angular size to 
be < 3.8 mas (or < 30 pc at z abs ). Given the compact na- 
ture of the radio source it is intriguing to see the difference 
between the H2 and 21-cm absorptions. This clearly shows 
that the absorbing gas contains a mixture of different phases 
at parsec scales. 



7 RESULTS AND DISCUSSION 

We have carried out a systematic search for 21-cm absorp- 
tion in 10 DLAs at z abs > 2 using GMRT and GBT. We de- 
tect 21-cm absorption in only one of them. From our sam- 



ple we find the 21-cm detection rate is 13% for a J rdv 
limit of 0.4 km/s (the detection limit reached in the case of 
J1337+3152). We also obtained 1420 MHz VLBI images for 
the sources in our sample. 

The 21-cm detection at z > 2 seems to favour sys- 
tems with high met a llicity and/or high N(H 1) (see also 
iKanekar et~ai1 |2009| ; ICurran et afl |2010| ). This basically 
means that the probability of detecting cold components 
that can produce detectable 21-cm absorption is higher in 
systems with high values of JV(H 1) and Z. However, we do 
not find any correlation between the integrated optical depth 
(or Ts// C ) and N(R 1) or metallicity. 

It is important to address the covering factor issue be- 
fore drawing any conclusions on Ts- Ideally one should do 
high spatial reso lution VLBA spe ctroscopy for this purpose 
(see for example lLane et al.ll200ol ). However, this is not pos- 
sible at present specially for z > 2 absorbers. Therefore, 
we proceed by assuming that the core fraction found in the 
VLBA images as the covering factor of the absorbing gas 
(as in the case of IKanekar et all 120091 ). We find that more 
than 50% of DLAs have weighted mean spin temperature 
(Ts) in excess of 700 K. For the assumed temperature of the 
CNM gas Tg = 200 K (as seen in H2 components in high-z 
DLAs) we find that more than 73% of H I in such systems 
is originating from WNM. The median value CNM fraction 
(i.e /(CNM)) obtained for the detections and the median 
value of upper limits in the case of non-detections are in the 
range 0.2 to 0.25. 

We study the connection between 21-cm and H2 absorp- 
tion in a sub-sample of 13 DLAs where both these species 
can be searched for. We report the detection and detailed 
analysis of H2 molecules in the z abs =3.3871 DLA system 
towards J0203+1134 where 21-cm absorption is also de- 
tected. For a b parameter in the range 1-5 km s _1 we find 
14.57<log Af(H2)<16.03. The inferred kinetic temperature 
is in the range 48-108 K based on Tbi of H2. However no 
21-cm absorption is detected at the very position of this H2 
component. This suggests that the H I column density as- 
sociated with this component is < 10 19 cm -2 . However, the 
lack of proper coincidence between 21-cm and any of the 
strong UV absorption components may also mean that the 
radio and optical sight lines probe different volumes of the 
gas. 

In the case of 8 DLAs, neither 21-cm nor H2 are de- 
tected. Typical upper limits on the molecular fraction (/h 2 ) 
in these systems are < 10 -6 . The lack of H 2 in DLAs can 
be explained if the H I gas originates from low density re- 
gions photoionized by the metagalactic UV (see for ex ample, 
iPetitiean et~al]| 19921 [20001 ; iHirashita fc Ferrarall2005l ). This 
also indi cates that the volume filli ng factor of H2 in DLAs 
is small (|Zwaan fc Prochaska|[2006l ). Typical limits obtained 
for Ts in these systems are consistent with only a small 
fraction of the H I gas originating from the CNM phase as 
suggested by the lack of H2 absorption. 

In two cases strong H2 absorption is detected and ki- 
netic temperatures are in the range 100-200 K, but 21- 
cm absorption is not detected. Even in two cases where 
both the species are detected they do not originate from 
the same velocity component. The lack of 21-cm absorp- 
tion directly associated with H2 indicates that only a small 
fraction (typically < 10%) of the neutral hydrogen seen in 
the DLA is associated with the H2 components (see also 
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iNoterdaeme et all l20ld ) . This implies that the molecular 
fractions /(H2) reported from the H2 surveys should be con- 
sidered as conservative lower limits for the H2 components. 

For two of the H2-bearing DLAs with density measure- 
ments based on C I fine-structure excitation we derive an 
upper limit on the line of sight thickness of < 15 pc. This 
is consistent with the size estimate for the Lb-bearing gas 
in z a bs= 2.2377 DLA towards Q1232+082 based on partial 
coverage (|Balashev et al.ll201ll ). 

In principle, the presence of H2 and 21-cm absorptions 
in a single component provides a unique combination to si- 
multaneously constrain the variation of the fine-structure 
constant (a), the electron-to-proton mass ratio (fi) and the 
proton G-factor. As shown here, DLAs with 21-cm and H2 
detections are rare. Even in these cases the presence of mul- 
tiphase structure at parsec scale is evident, introducing ve- 
locity shifts between the different absorption components 
that will affect the constraints on the variation of constants. 
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